arXiv: 1507.06319vl [astro-ph.SR] 22Jul2015 


Draft version July 24, 2015 

Preprint typeset using DTfrjX style emulateapj v. 5/2/11 


A SECOND CASE OF OUTBURSTS IN A PULSATING WHITE DWARF OBSERVED BY KEPLER 

J. J. Hermes 1 , M. H. Montgomery 2 , Keaton J. Bell 2 , P. Chote 1 , B. T. Gansicke 1 , Steven D. Kawaler 3 , 
J. C. Clemens 4 , B. H. Dunlap 4 , D. E. Winget 2 , and D. J. Armstrong 1 

Draft version July 24, 2015 

ABSTRACT 

We present observations of a new phenomenon in pulsating white dwarf stars: large-amplitude 
outbursts at timescales much longer than the pulsation periods. The cool ( T e g = 11,010 K), hydrogen- 
atmosphere pulsating white dwarf PG 1149+057 was observed nearly continuously for more than 78.8 d 
by the extended Kepler mission in K2 Campaign 1. The target showed 10 outburst events, recurring 
roughly every 8 d and lasting roughly 15 hr, with maximum flux excursions up to 45% in the Kepler 
bandpass. We demonstrate that the outbursts affect the pulsations and therefore must come from 
the white dwarf. Additionally, we argue that these events are not magnetic reconnection flares, and 
are most likely connected to the stellar pulsations and the relatively deep surface convection zone. 

PG 1149+057 is now the second cool pulsating white dwarf to show this outburst phenomenon, after 
the first variable white dwarf observed in the Kepler mission, KIC 4552982. Both stars have the same 
effective temperature, within the uncertainties, and are among the coolest known pulsating white 
dwarfs of typical mass. These outbursts provide fresh observational insight into the red edge of the 
DAV instability strip and the eventual cessation of pulsations in cool white dwarfs. 

Subject headings: stars: individual (PG 1149+057)-stars: white dwarfs-stars: oscillations (including 
pulsations)-stars: variables: general-stars: evolution-stars 


1. INTRODUCTION 

The vast majority of all stars in our Galaxy will end 
up as white dwarfs, which provide important observa¬ 
tional boundary conditions on the final stages of stellar 
evolution. White dwarfs have electron-degenerate cores 
overlaid by a thin, non-degenerate atmosphere, and pas¬ 
sively release their residual thermal energy. When a DA 
(hydrogen-atmosphere) white dwarf eventually cools be¬ 
low 12,500 K, it develops a hydrogen partial-ionization 
zone at its surface. This superficial convection zone 
throttles heat transport and efficiently drives pulsations 
(Brickhill 1991a; Goldreich & Wu 1999a). 

The flux changes in variable DA white dwarfs (DAVs, 
aka ZZ Cetis) are consistent with surface temperature 
variations caused by non-radial g-mode pulsations. By 
matching the observed pulsation periods to theoretical 
models, we can reveal considerable detail about the in¬ 
terior of white dwarfs, including their overall mass, tem¬ 
perature, the mass of their outer hydrogen and helium 
envelopes, their rotation, and even characteristics of their 
cores (see reviews by Winget & Kepler 2008; Fontaine & 
Brassard 2008; Althaus ct al. 2010). 

As white dwarfs cool their convection zones deepen. 
Since the oscillation modes excited are characterized by 
the thermal timescale of the driving zone, a deeper con¬ 
vection zone excites long-period pulsations, up to an 
observed maximum of about 1200 s for a 0.6 Mq white 
dwarf (Mukadam et al. 2006). Empirical studies of the 

j. j. hermes@ Warwick, ac. uk 

1 Department of Physics, University of Warwick, Coventry- 
CV4 7AL, UK 

2 Department of Astronomy, University of Texas at Austin, 
Austin, TX-78712, USA 

3 Department of Physics and Astronomy, Iowa State Univer¬ 
sity, Ames, IA-50011, USA 

4 Department of Physics and Astronomy, University of North 
Carolina, Chapel Hill, NC - 27599-3255, USA 


coolest members of the DAV instability strip are com¬ 
plicated by the relatively long-period and less coherent 
pulsations (Kanaan et al. 2002). However, the Kepler 
space telescope provides an unprecedented tool for long¬ 
term monitoring of the coolest pulsating white dwarfs. 

Ten pulsating white dwarfs were eventually discovered 
late in the original Kepler mission, but only two were 
observed for more than a month before the failure of 
the second reaction wheel in 2013 May (Hermes et al. 
2011; Greiss et al. 2014). Curiously, the first and longest- 
studied DAV observed by Kepler, KIC 4552982, showed 
a completely unprecedented behavior: large-amplitude 
flux excursions which raised the overall brightness by 
up to 17%, recurred on average every 2.7 d and lasted 
4 — 25hr in duration (Bell et al. 2015). These outbursts 
recurred on a much longer timescale than the pulsations. 

The extended Kepler mission ( K2 ) affords an excep¬ 
tional opportunity to extend space-based monitoring for 
up to 85 d to a plethora of new (and brighter) white 
dwarfs (Howell et al. 2014). Already K2 has illuminated 
frequency variability in the cool DAV GD 1212 (Hermes 
et al. 2014) and a hot DAV descended from common- 
envelope evolution (Hermes et al. 2015). 

The pulsating white dwarf PG 1149+057 was observed 
for more than 78.8 d in K2 Campaign 1, and we present 
here an analysis of the second case of sporadic, high- 
amplitude outbursts in a cool DAV. In this Letter we 
present our K2 observations, establish that the outbursts 
are happening on the white dwarf, and analyze the out¬ 
bursts and their relation to the stellar pulsations. 

2. K2 CAMPAIGN 1 TIME-SERIES PHOTOMETRY 

There is little high-speed photometry on PG 1149+057, 
despite the fact that there are just 15 brighter known 
DAVs. Pulsations were announced using one very short 
(< 0.7hr) run in 2005 May by Voss et al. (2006), who 
present the only previously published light curve. 
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Fig. 1.— Representative portions of the K2 Campaign 1 light curve of the pulsating white dwarf PG 1149+057. The top left panel shows 
the first 25 d of observations; three outburst events are denoted in green. The bottom left panel shows 7.2 hr of data on the second day of 
K2 observations; the white dwarf pulsations are clearly visible, and underplotted is a best-fit to the three highest-amplitude signals (with 
periods of 1145.7s, 998.1s, and 1052.8s). The right panel shows 7.2hr during the second outburst, with points connected in green. 


Our K2 Campaign 1 coverage of PG 1149+057 
(EPIC 201806008, K p = 15.0 mag) was requested in 
short-cadence (exposures every 58.8 s) and spans from 
2014 June 3 00:00:42 UT to 2014 August 20 20:16:36 
UT. Every 6 hr, the Kepler spacecraft checks its pointing 
and, if it needs correction, fires its thrusters. This can 
cause significant discontinuities in the light curve. We 
have extracted and detrended the photometry using the 
tools outlined in Armstrong et al. (2015), which compen¬ 
sates for artifacts in the light curve caused by attitude 
corrections. 

Our final data set has 107,443 points, yielding a duty 
cycle >92.8%. Figure 1 shows the first 25 d of K2 cov¬ 
erage, including the first three of a total of 10 out¬ 
bursts. Inspecting the light curve in more detail reveals 
the 1 — 3% peak-to-peak variability caused by the white 
dwarf pulsations. 

We display the entire K2 light curve in the top panel of 
Figure 2, marking each outburst in green. The multi-day 
gap beginning near 2456849 BJD was caused by a ma¬ 
jor pointing shift in the spacecraft, which subsequently 
added minor, long-timescale variations about the mean 
flux level. 

We have inspected movies of the Kepler images around 
each of these outbursts using the k2flix package (Bar- 
entsen 2015) and see no evidence for passing solar sys¬ 
tem bodies causing brightness increases. Kepler has a 
large plate scale (4 arcsec pixel” 1 ), but the nearest source 
with r < 19.5 mag is SDSS J115155.09+052959.3, which 
is more than 1.3arcmin to the north, well outside the 
extracted target pixels. 

3. ARCHIVAL SPECTROSCOPY AND PHOTOMETRY 


We describe here the time-averaged atmospheric pa¬ 
rameters for PG 1149+057, and exclude the possibility 
that this white dwarf has a flaring, low-mass companion. 

To date, the most consistent picture of the parame¬ 
ters of white dwarfs in and around the DAV instability 
strip has been undertaken by Gianninas et al. (2011). 
Their one-dimensional, ML2/a = 0.8 atmospheric pa¬ 
rameters place PG 1149+057 near the empirical red edge 
of the DAV instability strip, with T e g = 11,360 ± 170 K 
and log g = 8.21 ± 0.05. Tremblay ct al. (2013) intro¬ 
duced terms to correct for the three-dimensional (3D) 
dependence of convection, which refine the final deter¬ 
minations for PG 1149+057 to T eS = 11,060 ± WOK 
and log g = 8.06 ± 0.05. This corresponds to a mass 
of 0.64 ± 0.03 Mq and a distance of 39.3 ± 2.4pc using 
the mass-radius relationship of Renedo et al. (2010). 

High-resolution spectra taken as part of the ESO Su¬ 
pernovae Type la Progenitor Survey show no evidence 
of metals from ongoing accretion (Voss et al. 2006), and 
the radial velocities of this white dwarf do not vary to 
a limit of 4.0kms” 1 over 4hr (Maxted & Marsh 1999). 
We have inspected the Wide Angle Search for Planets 
survey light curve of PG 1149+057, which covers more 
than 3yr but is at the faint end for the survey. Despite 
binning the data into 30-min intervals, the la standard 
deviation of all points is 0.17mag, so we cannot signif¬ 
icantly constrain whether outbursts have been observed 
from previous ground-based studies. 

We have examined the wealth of archival photometry 
on PG 1149+057 and show in Figure 3 that the spectral 
energy distribution (SED) is consistent with a single, iso¬ 
lated white dwarf from the near-ultraviolet out to the 
mid-infrared. The SED uses photometry from GALEX- 
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Fig. 2. — The top panel shows the entire K2 light curve of PG 1149+057. We mark the 10 outburst events in green and show the data 
in quiescence as black points. The bottom panel shows FTs both in (green) and out of outburst (black). Both subsets have pulsations 
excited above the significance threshold, described in the text and shown as dark green dashed-dotted and blue dashed lines, respectively. 
Pulsations persist during outburst but have higher amplitudes and shorter periods than in quiescence. 


NUV (Morrissey et al. 2007), SDSS-ugri;? and UKIDSS- 
JHK (Girven et al. 2011), and Spitzer 3.6- and 4.5 pm 
photometry (Barber et al. 2012). Underplotted in Fig¬ 
ure 3 is an 11,000 K, log<? = 8.0 DA white dwarf model 
atmosphere computed by Koester (2010), normalized to 
the SDSS-r band. 

The lack of infrared excess as compared to our white 
dwarf model rules out any possible companion ear¬ 
lier than L4 using the data of Cushing et al. (2005). 
Our K2 light curve is therefore composed entirely of 
PG 1149+057, and the outbursts observed are the result 
of brightening events on the surface of the white dwarf. 



Fig. 3. — The spectral energy distribution from archival data 
demonstrates that there is no companion to PG 1149+057 respon¬ 
sible for the outbursts. The data are consistent with an isolated 
white dwarf (this an 11,000 K model) from the near-ultraviolet out 
to the mid-infrared (see text for sources). The 3<x photometric 
uncertainties are smaller than each point. 


4. DISCUSSION OF THE OUTBURSTS 
4.1. Outburst Characteristics 

The outburst events shown in the top panel of Figure 2 
confirm a new phenomenon observed in pulsating white 
dwarfs. Only one other pulsating white dwarf displays 
similar recurring behavior: The cool DAV KIC 4552982 
(WD J1916+3938), observed for nearly 2.5 yr in the orig¬ 
inal Kepler mission (Bell et al. 2015). 

That DAV has very similar 3D atmospheric parame¬ 
ters to PG 1149+057 - T e g = 10,860 ± 120 K and logg 
= 8.16 + 0.06 — and it also has a very similar pulsa¬ 
tion spectrum. The only other reminiscent event seen in 
a pulsating white dwarf occurred during monitoring of 
the cool variable DB (helium-atmosphere) white dwarf 
GD 358, which in 1996 underwent a dramatic change in 
excited frequencies accompanied by a rapid increase in 
fractional amplitude, its sforzando (Montgomery et al. 
2010 ). 

The outbursts in KIC 4552982 increase its mean flux 
by between 2 — 17%, last 4 — 25 hr, and recur on an av¬ 
erage every 2.7d (Bell et al. 2015). Everything about 
the outbursts in PG 1149+057 is scaled up by roughly a 
factor of three. The outbursts here increase the quies¬ 
cent flux by between 5 — 45%, have durations on aver¬ 
age of 14.5hr (ranging from 9.3 — 36.4hr), and recur on 
a timescale between 3.5 — 11.5 d, on average every 8d. 
These rough timescales were determined by eye, noting 
when the point-to-point scatter rose and fell to the qui¬ 
escent values. With so few events, an autocorrelation 
function reveals only that most events last more than 
12 hr. There does not appear to be any periodicity to 
the outbursts. 
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Fig. 4.— Diagnostics of the pulsations before (red), during (green) and after (blue) the seventh outburst observed on PG 1149+057. The 
top panel shows the K2 light curve, while the bottom panels show the FT of each subset. Signals exceeding the significance thresholds (see 
text) are marked with appropriate-colored dots and show differences in the excited pulsation spectra. Amplitudes are essentially all lower 
after outburst, but subsequently grow (see Figure 5). Significant pulsations are also visible during the brightening event (green squares), 
and appear at shorter periods (higher frequencies) and higher amplitudes. The top right panel shows the window function of each subset. 


4.2. Outbursts Affect the Pulsations 

We see convincing evidence that the pulsations persist 
on the white dwarf during the outburst, and actually 
change properties relative to those excited in quiescence. 
It is important to note that if these outbursts came from 
a contaminating source the pulsation amplitudes would 
decrease as a result of flux dilution. That the pulsation 
amplitudes actually increase proves that the outbursts 
are on the white dwarf. 

The bottom panel of Figure 2 shows Fourier transforms 
(FTs) of the entire light curve, differentiated by data in 
and out of outburst. Both datasets have peaks in the 
FT above the respective 3er significance threshold, which 
was determined by randomly shuffling the fluxes for each 
point in the light curve and computing the highest peak 
in the entire resultant FT for 10,000 random permuta¬ 
tions (see Hermes et al. 2015). 

Figure 2 shows that the pulsations in outburst have 
shorter periods and higher amplitudes than those in qui¬ 
escence; the unweighted mean period before, during, and 
after the seventh outburst goes from 1083.3s to 919.7s 
to 1072.7s, respectively. If the flux increase during the 
outbursts are due to an increase in the average effective 
temperature of the star, the formalism of Wu & Goldreich 
(1999) predicts that the surface amplitudes of the modes 
should increase, with the enhancement being greater for 
shorter-period modes. The 1000 s inodes are expected to 
have growth times of order days (Goldreich & Wu 1999a). 

Smoothing the light curve by 3000 s, over several pul¬ 
sation cycles, the outbursts increase the mean flux by 
up to 14%, which would require a global 750 K temper¬ 
ature increase from the quiescent effective temperature 


of 11,060 K. The pulsations induce even more dramatic, 
localized temperature excursions. 

Figure 4 explores the light curve surrounding and dur¬ 
ing the seventh observed outburst, between Days 46 — 56 
from the start of K2 monitoring, using an equal amount 
of data before and after the outburst. An FT of each 
subset clearly shows that pulsations are excited to sig¬ 
nificant amplitudes, even during the outburst. 

Many of the significant signals in outburst are not 
present in quiescence, especially the mode at 949.3 ± 
0.8/iHz (1053.4 ± 0.8s), which decreases in amplitude 
the longer we monitor it after outburst. Many modes 
appear to eventually return at roughly the same period 
after the outburst. However, the amplitude evolution of 
these pulsations holds key insight. 

In essentially all cases in Figure 4, the amplitudes 
shown in the subset immediately after the outburst are 
lower than the pre-outburst subset. There are more data 
in this post-outburst subset, and we can clearly see many 
modes growing after the outburst. We explore this am¬ 
plitude evolution using a running FT, shown in Figure 5. 
We have computed this diagram by taking an FT across a 
sliding 3-day window, stepped every hour with 0.05 /rHz 
frequency resolution. 

This time-dependent FT demonstrates that the out¬ 
burst has depleted considerable pulsation energy. The 
apparent resetting of amplitudes offers the unprece¬ 
dented opportunity to actually measure the linear growth 
rates of pulsation modes in a white dwarf, although that 
is outside the scope of this Letter. 

Finally, the bottom right panel of Figure 4 shows the 
relatively stable group of three symmetric signals cen¬ 
tered at 2425.0/zHz (the m = ±1 components are visible 
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Fig. 5. — Running FT of the K2 light curve near the seventh outburst shown in Figure 4. Vertical dotted lines mark the maximum 
brightness during outburst; our 3-day sliding window smears the events. This plot tracks the post-outburst frequencies, many of which 
have growing amplitudes until the onset of the eighth outburst peaking near Day 61. 


in this short subset but not the central component). This 
pulsation triplet arises from rotation of the white dwarf 
(Pesnell 1985), from which we measure an £ = 1 fre¬ 
quency splitting of roughly 4.7 /iHz, which corresponds 
to a rotation period of roughly 1.2 d in the asymptotic 
limit. We reserve a full period determination and de¬ 
tailed asteroseismic analysis of the many pulsations in 
PG 1149+057 for a future publication. 

4.3. Consideration of Possible Outburst Mechanisms 

We consider various timescales in white dwarfs to ex¬ 
plore different underlying mechanisms for the outbursts. 

White dwarfs have very short dynamical timescales, 
roughly a few seconds. While this is a global estimate, we 
expect all processes which are not in hydrostatic equilib¬ 
rium to locally occur on this or an even shorter timescale. 
The timescales for the duration and recurrence of the out¬ 
bursts are clearly much longer than this. This likely rules 
out magnetic reconnection events as responsible for the 
outbursts, since the rise times of magnetic flares occur 
at roughly the dynamical timescale, and the rise times of 
the outbursts here are of order hours. 

On the other hand, thermal processes can occur on 
longer timescales. To explore this more quantitatively, 
we have computed a theoretical white dwarf envelope 
appropriate to this star (11,000K, logg = 8.0, ML2/a = 
1.0) using the Warsaw - New Jersey stellar envelope code 
(Pamyatnykh 1999). We find that the thermal timescale 
corresponding to the average duration of the outbursts, 
14.5 hr, occurs at a depth of 10 -11 ' 5 of the star by mass, 
placing it slightly deeper but near the base of the con¬ 
vection zone. 

We note that the recurrence timescale of roughly 8 d is 
mapped to a deeper layer, the outer 10“ 10 of the star by 
mass. It is possible that processes at these depths play a 
role in the outbursts. However, the recurrence timescale 
probably rules out runaway thermonuclear events as the 
source of the temporary brightness increases, since the 
equilibrium temperature in this region is below 10 6 K, 


significantly too low for nuclear burning. 

An occasional deposition of asteroids onto the surface 
of the white dwarf could contribute the energy required 
to match the brightness increases, on average roughly 
10 34 erg integrated over each event, and we now have 
compelling evidence that disrupted planetesimals fre¬ 
quently pollute the atmospheres of white dwarfs (e.g., 
Koester et al. 2014). However, this pollution causes 
metal absorption, which we can rule out from spec¬ 
troscopy (Voss et al. 2006). 

Since the only two known outbursting white dwarfs 
have very similar effective temperatures and pulsation 
spectra, it is logical to suspect a connection between the 
events and the pulsations and the deep convection zone 
that drives them, especially since we have demonstrated 
here that the pulsations change before, during, and after 
outbursts. 

We consider a model where the outbursts represent 
a temporary, rapid reassignment of kinetic energy away 
from pulsations. It is possible to transfer energy be¬ 
tween pulsation modes through resonant mode coupling 
(Dziembowski 1982), which has previously been invoked 
to explain the amplitude evolution of 5 Scuti stars in the 
Kepler mission (Breger & Montgomery 2014). We note 
that there is still considerable kinetic energy in the modes 
which are damped in radiative regions of the star as well 
as from turbulent dissipation in the convection zone. 

A theoretical framework exists for parametric instabil¬ 
ity via mode coupling of white dwarf pulsations, in which 
energy is rapidly channelled from an observed, overstable 
parent mode to damped daughter modes (Wu & Goldre- 
ich 1999). It is possible that the energy of a mode (or 
multiple modes) grows linearly until it reaches a critical 
threshold, after which it enters a nonlinear regime and 
rapidly transfers its energy to resonant daughter modes, 
or a cascade of resonant daughter modes. These daugh¬ 
ter modes may be quickly damped by turbulence in the 
convection zone and thus deposit their newfound energy 
there. However, more detailed energetics calculations are 
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required to determine if there is sufficient energy in one 
(or even multiple) modes to power the observed luminos¬ 
ity increases during outburst. 

4.4. Conclusions 

We present confirmation of a new phenomenon in pul¬ 
sating white dwarfs by analyzing the second case of out¬ 
bursts in a cool DAV observed by the Kepler spacecraft, 
PG 1149+057. We suggest that these outbursts are con¬ 
nected with the pulsations, perhaps as a result of non¬ 
linear mode coupling. 

Still, exciting questions remain. Why have we not pre¬ 
viously seen outbursts in extensively studied cool DAVs, 
such as GD 154 and G29-38 (Pfeiffer et al. 1996; Klein- 
man et al. 1998)? How do the horizontal surface veloci¬ 
ties change from quiescence to outburst (e.g., van Kerk- 
wijk et al. 2000)? If resonant mode coupling of pulsations 
is responsible, what mechanism triggers the outburst? 
Are the outbursts responsible for the eventual cessation 
of pulsations and the empirical red edge of the DAV in¬ 
stability strip? 

Theory can reasonably reproduce the onset of pulsa¬ 
tions observed at the blue edge of the DAV instability 
strip (Brickhill 1991b; Wu & Goldreich 1999). However, 
there is still not a clear picture for the red edge where 
pulsations shut down, which is observed empirically to 
occur at roughly 10,500 K for a canonical 0.6 Mq white 
dwarf (Gianninas et al. 2014). Current non-adiabatic cal¬ 
culations predict a much cooler red edge, below 6000 K 
(Van Grootel et al. 2012). 

Hansen et al. (1985) proposed that longer-period 
modes would not be reflected off the surface, setting 
a condition for the red edge. Subsequently, Brickhill 
(1991b) and Goldreich & Wu (1999b) suggested that tur¬ 


bulent viscosity in the convection zone may damp the 
longest-period modes and thus define the red edge. How¬ 
ever, neither theory provides a complete picture. Per¬ 
haps this newly discovered outburst phenomena plays 
some role in shutting down pulsations in cool DAVs. 

Long-timescale monitoring is necessary to catch and 
constrain future outbursts from cool pulsating white 
dwarfs, and it is a testament to the utility of the Kepler 
instrument to non-planetary science that these events 
were first discovered thanks to its unblinking gaze. For¬ 
tunately, the extended Kepler mission continues its tour 
of new K2 fields along the ecliptic every three months, 
offering the possibility to survey additional cool DAVs 
in the coming years. We look forward to future theo¬ 
retical and observational constraints to this unexpected 
phenomenon. 


We acknowledge useful discussions with Tom Marsh, 
Jim Fuller, and Yanqin Wu, and a thoughtful report 
from an anonymous referee. We thank Simon Walker 
for sharing the WASP light curve of PG 1149+057 and 
thank Detlev Koester for making his DA atmospheric 
models freely available. J.J.H., P.C. and B.T.G. acknowl¬ 
edge funding from the European Research Council under 
the European Union’s Seventh Framework Programme 
(FP/2007-2013) / ERC Grant Agreement n. 320964 
(WDTracer). K.J.B., M.H.M. and D.E.W. gratefully ac¬ 
knowledge the support of the NSF under grant AST- 
1312983. M.H.M. acknowledges the support of NASA 
under grant NNX12AC96G. Funding for the Kepler mis¬ 
sion is provided by the NASA Science Mission Direc¬ 
torate. 


REFERENCES 


Althaus, L. G., Corsico, A. H., Isern, J., & Garcfa-Berro, E. 2010, 
A&A Rev., 18, 471 

Armstrong, D. J., Kirk, J., Lam, K. W. F., et al. 2015, A&A, 579, 
A19 

Barber, S. D., Patterson, A. J., Kilic, M., et al. 2012, ApJ, 760, 26 
Barentsen, G. 2015, Astrophysics Source Code Library, 1503.001 
Bell, K. J., Hermes, J. J., Bischoff-Kim, A., et al. 2015, 
arXiv: 1506.07878 

Breger, M., & Montgomery, M. H. 2014, ApJ, 783, 89 
Brickhill, A. J. 1991a, MNRAS, 251, 673 
Brickhill, A. J. 1991b, MNRAS, 252, 334 

Cushing, M. C., Rayner, J. T., & Vacca, W. D. 2005, ApJ, 623, 
1115 

Dziembowski, W. 1982, AcA, 32, 147 
Fontaine, G., & Brassard, P. 2008, PASP, 120, 1043 
Gianninas, A., Bergeron, P., & Ruiz, M. T. 2011, ApJ, 743, 138 
Gianninas, A., Dufour, P., Kilic, M., et al. 2014, ApJ, 794, 35 
Girven, J., Gansicke, B. T., Steeghs, D., & Koester, D. 2011, 
MNRAS, 417, 1210 

Goldreich, P., & Wu, Y. 1999a, ApJ, 511, 904 
Goldreich, P., & Wu, Y. 1999b, ApJ, 523, 805 
Greiss, S., Gansicke, B. T., Hermes, J. J., et al. 2014, MNRAS, 
438, 3086 

Hansen, C. J., Winget, D. E., & Kawaler, S. D. 1985, ApJ, 297, 
544 

Hermes, J. J., Mullally, F., 0stensen, R. H., et al. 2011, ApJ, 741, 
L16 

Hermes, J. J., Charpinet, S., Barclay, T., et al. 2014, ApJ, 789, 85 
Hermes, J. J., Gansicke, B. T., Bischoff-Kim, A., et al. 2015, 
MNRAS, 451, 1701 

Howell, S. B., Sobeck, C., Haas, M., et al. 2014, PASP, 126, 398 


Kanaan, A., Kepler, S. O., & Winget, D. E. 2002, A&A, 389, 896 
Kleinman, S. J., Nather, R. E., Winget, D. E., et al. 1998, ApJ, 
495, 424 

Koester, D. 2010, In Memorie della Societa Astronomica Italiana, 
81, 921 

Koester, D., Gansicke, B. T., & Farihi, J. 2014, A&A, 566, A34 
Maxted, P. F. L., & Marsh, T. R. 1999, MNRAS, 307, 122 
Montgomery, M. H.. Provencal, J. L., Kanaan, A., et al. 2010, 
ApJ, 716, 84 

Morrissey, P., Conrow, T., Barlow, T. A., et al. 2007, ApJS, 173, 
682 

Mukadam, A. S., Montgomery, M. H., Winget, D. E., Kepler, 

S. O., & Clemens, J. C. 2006, ApJ, 640, 956 
Pamyatnykh, A. A. 1999, AcA, 49, 119 
Pesnell, W. D. 1985, ApJ, 292, 238 

Pfeiffer, B., Vauclair, G., Dolez, N., et al. 1996, A&A, 314, 182 
Renedo, I., Althaus, L. G., Miller Bertolami, M. M., et al. 2010, 
ApJ, 717, 183 

Tremblay, P.-E., Ludwig, H.-G., Steffen, M., Sz Freytag, B. 2013, 
A&A, 559, AA104 

Van Grootel, V., Dupret, M.-A., Fontaine, G., et al. 2012, A&A, 
539, A87 

van Kerkwijk, M. H., Clemens, J. C., & Wu, Y. 2000, MNRAS, 
314, 209 

Voss, B., Koester, D., 0stensen, R., et al. 2006, A&A, 450, 1061 
Winget, D. E., & Kepler, S. O. 2008, ARA&A, 46, 157 
Wu, Y., & Goldreich, P. 1999, ApJ, 519, 783 
Wu, Y., & Goldreich, P. 2001, ApJ, 546, 469 


